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FILM  ELASTIC  PROPERTIES 
DETERMINED  BY  THE  INDENTATION 
TEST-THEORETICAL  CONSIDERATIONS 


INTRODUCTION 

Intense  interest  in  thin-film  technology  has  been  spurred  by  the  growing  importance  of  micro¬ 
electronics,  metal  matrix  composites,  ceramic  matrix  composites,  high  transition  temperature  super¬ 
conducting  films,  and  compositionally  modulated  materials.  Correspondingly,  the  determination  of 
film  proper  lies  has  become  of  greater  importance.  The  need  for  techniques  to  study  the  mechanical 
properties  of  thin  films,  which  are  usually  defined  as  coatings  of  thickness  of  up  to  a  few  microme¬ 
ters,  has  recently  rekindled  an  interest  in  nricrohardness  and  submicroindentation  devices.  The  con¬ 
tinuous  indentation  test,  also  known  as  the  impression  test,  measures  the  elastic  constants  of  bulk 
materials  [1-3].  In  these  tests,  the  displacement  and  the  load  of  the  indenter,  which  is  pressed  onto  a 
specimen  at  a  predetermined  speed,  is  recorded.  The  slope  of  the  unloading  curves  in  the  load  vs 
depth  plot  determines  the  elastic  modulus  of  the  specimen.  Recently,  a  high-resolution  nanoindenter 
was  used  to  determine  the  Young’s  modulus  of  thin  films  [4]  from  the  linear  unloading  portion  of  the 
indentation  load  vs  depth  curves.  It  was  also  shown  that,  as  the  indentation  depth  d  approaches  the 
dimension  of  the  film  thickness  h ,  the  influence  of  the  substrate  can  be  detected  because  of  the  chang¬ 
ing  contributions  of  the  film  and  the  substrate  to  the  measured  elastic  constant.  For  small  indentation 
depth  relative  to  the  film  thickness,  the  data  approach  that  expected  for  bulk  film  material  since  the 
indenter  interrogates  the  film  near  its  surface  only.  For  deep  indentations,  the  data  approach  that 
expected  for  bulk  substrate  material. 

The  indentation  problem  is  a  mixed  boundary  value  problem,  more  commonly  known  as 
Boussinesq’s  problem.  Harding  and  Sneddon  [S],  Sneddon  [6],  and  Miki  [7]  have  considered  this 
problem  when  a  semi-infinite  space  is  indented  by  a  cone,  a  sphere,  and  a  flat-ended  cylindrical 
punch.  The  flat-ended  cylindrical  punch  problem  for  an  elastic  layer  resting  frictionlessly  on  a  rigid 
foundation  has  been  considered  by  Lebedev  and  Ufliand  [8].  The  rigid  substrate  supporting  the  elas¬ 
tic  layer  is,  of  course,  a  mathematical  artifice  that  simplifies  the  analysis.  Some  investigators  con¬ 
sidered  the  indentation  of  an  elastic  layer  perfectly  bonded  to  an  elastic  half-space  made  up  of  dif¬ 
ferent  materials.  Wu  and  Chiu  [9]  considered  the  plane  strain  problem  and  reduced  the  mixed  bound¬ 
ary  value  problem  to  a  single  Fredholm  integral  equation  of  the  second  kind.  Dhaliwal  [10]  was  able 
to  reduce  the  axisynunetric  problem  of  a  flat-ended  cylindrical  punch  on  a  layered  elastic  medium  to  a 
Fredholm  integral  equation,  which  he  solved  approximately  to  estimate  the  safety  of  foundations  sup¬ 
porting  cylindrical  columns.  Subsequently  Dhaliwal  and  Rau  [11]  extended  the  analysis  of  Dhaliwal 
to  punches  of  arbitrary  profile  but  did  not  present  any  numerical  results. 

A  second  approach  to  the  solution  of  mixed  boundary  value  problems  is  to  replace  the  exact 
boundary  condition  by  an  approximate  one  and  solve  the  new  elasticity  problem.  Chen  and  Engel 
[12]  used  this  approach  to  study  the  impact  and  contact  stresses  caused  by  flat-ended  cylindrical  and 
parabolic  punches  on  a  composite  medium  consisting  of  one  or  two  elastic  layers  perfectly  bonded  to 
each  other  and  to  an  elastic  homogeneous  half-space. 
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The  punch  problem  (for  any  shape  punch)  for  an  elastic  layer  resting  frictionlessly  on  an  elastic 
substrate  (half-space)  has  not  been  considered.  The  punch  problems  for  conical  or  spherical  punches 
for  a  layer  perfectly  bonded  to  an  elastic  substrate  have  not  been  considered  either.  The  need  for  the 
solution  of  these  problems  for  determining  the  elastic  properties  of  thin  films  by  the  indentation  test  is 
obvious.  The  contact  between  the  thin  film  and  the  substrate  is  neither  frictionless  nor  perfectly 
bonded.  The  solutions  and  numerical  values  for  both  ideal  cases  should  be  obtained  to  set  guidelines 
for  the  indentation  test,  i.e.,  to  determine  what  the  proper  substrate  is  and  what  the  required  film 
thickness  is  so  that  the  film’s  elastic  constants  may  be  measured  within  a  predetermined  degree  of 
accuracy.  These  are  goals  of  this  investigation. 

FORMULATION  OF  THE  PROBLEM 

Consider  an  infinite  layer  of  thickness  h  with  elastic  constants  ni,  Pi,  overlaying  a  half-space 
substrate  with  elastic  constants  /X2>  >'2>  ^  axisymmetric  indenter  (punch)  with  a  flat,  conical, 

or  hemispherical  end  (Figs.  1  to  3,  respectively)  and  axis  normal  to  the  layer  surface,  pressing  fnc- 
tionlessly  on  the  layer.  Choose  cylindrical  coordinate  axes  (r,  ff,  z)  such  that  z  is  parallel  to  the  gen¬ 
eratrix  of  the  indenter,  r  is  perpendicular  to  z ,  d  is  the  angular  distance  between  a  reference  line  and 
r,  and  the  origin  of  coordinates  is  located  at  the  first  point  of  contact  between  the  indenter  and  the 
layer  (center  of  first  contact  area  for  the  flat-ended  indenter). 


t 

/ 


Fig.  1  —  A  composite  medium  consisting  of  an  elastic  layer  either  perfectly  bonded  to 
or  smoothly  overlaying  an  elastic  semi-infinite  substrate  indented  by  a  rigid  flat-ended 
cylindrical  indenter 
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P 


Fig.  2  —  A  composite  medium  consisting  of  an  elastic  layer  either  perfectly  bonded  to 
or  smoothly  overlaying  an  elastic  semi-infinite  substrate  indented  by  a  rigid  conical 
indenter 


Fig.  3  —  A  composite  medium  consisting  of  an  elastic  layer  either  perfectly  bonded  to 
or  smoothly  overlaying  an  elastic  semi-infinite  substrate  indented  by  a  rigid  spherical 
indenter 
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The  Papkovich-Neuber  functions  may  then  be  used  to  formulate  the  problem.  Thus,  in  terms  of 
a  pair  of  harmonic  functions  <fij  (r,  z)  and  'ij  (r,  z)  (/  =  1,  2),  the  displacement  and  stress  com¬ 
ponents  of  interest  may  be  written  as 

~  ^j.Z  ~ 


^j,zz  (^) 

"^rzi  ~  '^j.zz  ~  ^^j,rz' 


j  =  1,2,  kj  =  3  -  4Pj, 

where  the  letter  subscripts  following  a  comma  denote  differentiation  with  respect  to  the  indicated 
cylindrical  coordinates,  e.g.,  <pj^^  =  d^<Pj/drdz. 


The  harmonic  functions  <Pj{r,  z),  4fj(,r,  z)  may  be  expressed  in  terms  of  a  new  set  of  unknown 
functions  y4,  (X),  (/  =  1,  2,  . . . ,  6)  as  the  following  Hankel  integrals 

, »  JoC^r) 

(Aj  chXz  +  A2  sh\z)-^^-^  dX, 


Jo(Xr) 


^i('‘-z)  =  lo  (.43  sh  Xz  +  ^4  ch  Xz)  — —  (iX, 


(2) 


i'2(r,  z)  =  j"  Ase-^(^  -  ''%(Xr)  dX, 


oo 

>P2if,  z)  =  A(,e 


-Mz  -  A) 


Jo  (Xr) 


dX, 


where  ch  Xz  and  sh  Xz  are  hyperbolic  cosine  of  Xz  and  hyperbolic  sine  of  Xz  respectively  and  Jq  (Xr) 
is  the  Bessel  function  of  the  first  kind  of  order  0.  It  may  be  readily  be  seen  that  Eqs.  (2)  satisfy  the 
boundary  conditions  that  the  stresses  and  their  derivatives  in  both  media  vanish  for  r  >  0,  z  —00, 
and  for  z  >  0,  r  —  00 . 

Regardless  of  type  of  contact  between  layer  and  substrate  the  displacement  and  traction  bound¬ 
ary  conditions  of  the  surface  may  be  expressed  as 

«^,(r,  0)  =  d  -  S{r/a)  (0  <  r  <  d),  (3.1) 

<7ji(r,0)  =0  (a  <  r  <  00),  and  (3.2) 


Trjl(r,0)  =  0 


(0  <  r  <  00), 


(3.3) 
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where  d  is  the  depth  of  penetration  of  the  indenter,  a  is  the  radius  of  the  circle  of  contact,  and  the 
function  b{r /a)  is  prescribed  by  the  fact  that,  in  reference  to  the  tip  of  the  indenter  as  origin,  the 
punch  has  equation  z  =  b(r/a)  so  that  5(0)  =  0.  The  function  8(r/a)  for  conical,  hemispherical, 
and  flat-ended  cylindrical  indenters  on  the  half-space  have  been  given  by  Sneddon  [6]. 


With  the  functions  <Pj{r,  z)  and  ^jir,  z)  (j  =  1,2)  as  shown  in  Eq.  (2),  the  boundary  condi¬ 
tions  Eq.  (3.1)  to  (3.3)  may  be  expressed  by  using  Eq.  (1)  as  follows: 


,  00  /4i(X)  nx 

Jo  71^  W  = 


sh  Xh 


1  - 


[d  -  S(r/a)] 


,  »  A,(X)  -  2(1  -  .,)^2(X) . 

Jo  - Xn - >^o(Xr)  d\  =  0 


and 


Defining 


sh  Xh 


A^xCX)  =  (1  -  2vi)Ai  (X) 


(0  <  r  <  a). 


(a  <  r  <  oo), 


(0  <  r  <  oo). 


M(X)  = 


^i(X) 
sh  w 


,l  -  g(X)  = 


Ax(X) 


A^(X)  -  2(1  -  >',)/l2(X)  ’ 


(4.1) 


(4.2) 


(4.3) 


(5) 


where  w  =  Xh,  and  the  conditions  Eqs.  (4.1)  and  (4.2)  lead  to  a  system  of  two  integral  equations  for 
M(X): 

\^M(X)Jo(Xr)dX=f(r)  (0<r<a),  (6.1) 

and 

•'o(Xr)  dX  =  0  (a  <  r  <  oo),  (6.2) 

•’0  1  -  g  (X) 

where 


fir)  =  Id  -  dir /a)]. 

1  -  Vi 


(7) 


The  function  g(X)  can  be  obtained,  as  shown  in  detail  in  the  Appendix  A,  by  expressing  the 
unknown  functions  /1,(X)  (i  =  1,2,  ...,6)  in  terms  of  /4i(X)  for  different  boundary  conditions 
between  the  layer  and  the  semi-infinite  half-space.  This  is,  for  the  case  of  perfect  bonding,  continuity 
of  components  of  displacements  and  tractions  at  z  =  h  requires  that 


Urx  ir,h)  =  Ur2  (r,  h), 

ir,  h)  =  m,2  (r,  h),  (8) 

<7j,  (r,  5)  =  a,2  ir,  h),  and 

rrz\  ir,  h)  =  t^2  ir,  h). 
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for  0  <  r  s  00,  while  for  the  case  when  the  two  surfaces  are  in  smooth  contact  (frictionless),  con¬ 
tinuity  of  normal  components  of  displacement  and  tractions  and  the  conditions  of  vanishing  shear 
stress  components  at  the  interface  may  be  expressed  by 

(r,  h)  =  u^2  ('■.  h), 

a^i(r,h)  =  a^2  y  h),  md  (9) 

l-nl  (fy  *)  =  ^rz2  (fy  =  0, 

for  0  <  r  <  00 . 

The  solution  of  Eq.  (6)  is  sought  in  the  form  [9] 


M(X)  =  [1  -  g(X)l  cos  X/  dt. 


(10) 


where  <^(/)  is  the  solution  of  the  following  Fredholm  integral  equation  of  the  second  kind  with  a  con¬ 
tinuous  symmetrical  kernel: 

-  —  L  [G(5  +  /)  -t-  G(s  -  t)]<i>(s)ds  =  Fit)  (0  <  a),  (11) 

IT  O 

where 

G(jc)  =  g (X)  cos  Xx  rfX,  (12) 

=  -  1/(0)  +  t  l^f'it  sin  6)  dd],  (13) 

TT  U 

and 

fix)  =  ^/(^). 
dx 

After  determining  Fit),  the  unknown  function  </>(r)  can  be  solved  by  using  Eq.  (11).  Then  by 
using  Eq.  (10)  and  the  solutions  of  the  simultaneous  equations  of  i4,  (X),  Eqs.  (1)  and  (2)  give  a  com¬ 
plete  solution  to  the  contact  problem  under  consideration.  For  instance,  by  using  Eq.  (2)  and  the 
known  formula 


e  * 

JoiXr)  sin  \t  d\  =  0 


(0  <  t  <  r), 


(14) 


=  (f2  -  r2)-'/2 


it  >  r), 

it  is  easy  to  obtain  the  equation  for  the  distribution  of  the  normal  stresses  under  the  punch 


a,|(r,  0)  = 


0'(O 


(/"  -  n' 


dt  - 


<t>ia) 


ia^  -  ry 


(r  <  a). 


(15) 
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Substituting  Eq.  (7)  into  Eq.  (13)  we  have 


where 


F{t)  = 


»(1  -  I'll 


(16) 


Foir)  =  1 

and  7  =  t/a.  Introducing  the  dimensionless  quantities 

r  2a\d 

—  =  y  and  <^(r)  =  — - - -  H(j), 

a  t(1  -  V,) 

Eq.  (11)  then  assumes  the  form 

H(j) - ^  f  [A'(y  +  t)  +  K(y  -  r)]  H(y)  dy  =  FqW  (0  <  t  <  1), 

T  "'0 

where 


K{u)  =  X  Jn 
h  ■'0 


auw 

h 


dw. 


(17) 


(18) 


(19) 


(20) 


The  magnitude  of  the  applied  load  p  can  be  obtained  by  integrating  the  pressure  of  the  indenter 
on  the  layer  over  the  area  of  contact,  i.e.,  integrating  Eq.  (15)  over  the  area  of  the  circle  of  radius  a . 
Thus, 


P  = 


Aniod 
1  -  "1 


f  V(T)  dr. 


where  H{t)  is  known  after  Fq(t)  in  Eq.  (19)  is  determined. 


(21) 


The  function  Fgfr)  can  be  obtained  from  Eqs.  (9),  (13),  (16),  and  (17)  for  an  indenter  of  arbi¬ 
trary  axisymmetric  profile.  Intuitively,  the  contact  pressure  for  the  top  surface  of  the  layered  half¬ 
space  and  that  for  the  homogeneous  half-space  are  similar  in  many  respects  if  the  indenter  profile 
remains  the  same  in  both  cases.  For  example,  if  the  inder»er  has  a  sharp  comer,  the  stress  will  be 
singular  there.  If  the  identer  profile  is  smooth,  then  the  normal  stress  i  must,  on  physical  grounds, 
remain  finite  around  the  circle  r  =  a  (t  =  1).  This  continuity  of  normal  stresses  gives  an  additional 
equation 


=  H{\)  =  0, 


(22) 


which  is  a  consequence  of  Eqs.  (15)  and  (18).  Therefore,  in  the  case  of  indenters  with  smooth  pro¬ 
file,  such  as  conical  or  hemispherical,  Eqs.  (19)  and  (22)  have  to  be  solved  first  to  find  the  relation¬ 
ships  between  the  penetration  depth  d  and  the  radius  of  contact  area  a .  For  convenience,  let  us 
define  a  dimensionless  parameter  7,  which  is  a  function  of  h ,  ;ti,  V],  ^2,  >'2,  and  the  indenter  profile. 


7  = 


a 


an  ' 


(23) 
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where  0^  is  the  radius  of  the  contact  area  for  a  homogeneous  half-space  with  elastic  constants 
Hi  and  vi-  The  relationships  between  and  d  for  different  indenters  are  [5]; 

an  =  a,  (24) 


for  a  flat-ended  cylindrical  indenter  of  radius  a  (note  that  this  equation  is  independent  of  d); 


On  =  —  d  tan  a, 

T 

for  a  conical  indenter  with  included  angle  2a;  and 


Oh  =  2d 


In 


l±p 


(p  =  an/R), 


(25) 


=  (/?d)'''2  (a„  «  R),  (26) 

for  a  hemispherical  tip  indenter  of  radius  R . 

Since  when  using  the  indentation  test  to  determine  elastic  constants,  such  as  Young’s  modulus, 
shear  modulus,  or  Poisson  ratio,  each  constant  is  not  measured  directly  but  in  any  combinations  of 
pairs,  it  is  convenient  to  define  a  new  elastic  constant,  say 


f  = 


(27) 


which  may  be  called  the  impression  elastic  modulus,  or  the  impression  modulus  for  short.  Now  the 
relationship  between  p  and  z  for  each  type  of  indenter  may  be  obtained  as  follows. 


Flat-ended  Cylindrical  Indenter 

For  the  flat-ended  cylindrical  indenter  of  radius  a  (Fig.  1),  the  boundary  condition  Eq.  (3.1)  is 

Mj,  (r,  0)  =  d  (0  <  r  <  d). 


which  gives 

6(r/a)  =  0, /(r)  =  (0  <  r<a), 

and 

Fo(t)  =  1  (0  <  T  <  1),  (28) 

where  is  the  impression  modulus  for  the  homogeneous  half-space  with  elastic  constants  H\  and  v^. 
According  to  Eq.  (21),  the  relationship  between  the  p  and  the  d  can  be  written  as 
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r  =  ^xod  //(t)  dr,  (29) 

and 

Pi  =  tiod  (/=1.2).  (30) 

where  i  =  1  when  /i  —  oo  and  i  =  2  and  h  —  0. 

Conical  Indenter 

For  normal  penetration  by  a  rigid  cone  of  included  angle  2a  (Fig.  2),  the  boundary  condition, 
Eki-  (3.1),  gives 


6(r/a)  =  r  cot  a  (0  <  r  <a), 

then  Eqs.  (7)  and  (13)  become 

fir)  =  \^xd  1 
and 

F’o(r)  =  1  -  7(.T  (0  <  T  <  1), 

where  =  a  la^.  The  relationships  between  p ,  d,  and  a  are 

2  1  f ' 

p  =  —  7cfi^  tan  a  f  //  (r)  dr, 

TT  0 

=  y  cot  oc\^H  ir)dT, 

and 

Pi  =  tan  a  (i  =  1,  2). 


27c 


TTfl 


(0  <  r  <  a). 


(31) 


(32) 


(33) 


(34) 


Spherical  Indenter 


For  normal  penetration  by  a  rigid  spherical  indenter  of  radius  R  (Fig.  3),  the  same  as  before, 
the  strained  surface  of  the  elastic  layer  conforms  to  the  sphere  between  the  first  point  of  contact  and 
the  section  of  radius  a .  From  the  geometry  of  Fig.  3  we  find  that 


r 

-2 

Ml' 

=  R 

1  - 

1  - 

(0  <  r  <  a). 


(35) 


which  gives 


fir)  = 


1  - 


1  - 


1  - 


R- 


1/2 


7/p  In  (1  +  p)  -  In  (1  -  p) 


(0  <  r  ■<  a). 


(36) 
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and 


Fq{t)  =  1 


T  In  (1  +  pr)  —  In  (1  —  pr) 
7i  In  (1  +  p)  —  In  (1  -  p) 


(0  <  r  <  1), 


(37) 


where  7^  =  a /an  and  p  =  a/R.  If  the  indentation  is  small  (a  «  R  or  p  «  1),  then  approxi¬ 
mately 


and 


a  =  jsiRd)^^^  (a  «  R), 


(38) 


Fo(r)  =  1  -  7.V  (0  <  T  <  1,  p  «  1),  (39) 

and  from  Eqs.  (38)  and  (21) 

=  (40) 


P,  =  (1=1,2).  (41) 

The  formulation  of  the  problem  is  now  complete  for  each  of  the  three  types  of  indenter  considered. 
Numerical  results  for  some  cases  of  interest  are  presented  and  discussed  in  the  following  section. 

NUMERICAL  RESULTS 

The  mixed  boundary  value  problem  is  represented  by  the  Fredholm  integral  equation  of  the 
second  kind  as  shown  in  Eq.  (19).  For  a  given  function  k(u)  (Eq.  (20))  and  Fq{t)  (Eqs.  (28),  (31), 
or  (35)  depending  on  the  shape  of  the  indenter),  the  function  H{t)  is  solved  numerically  in  the  form 
of  a  Chebyshev  series  of  N  terms  (N  >  5)  [13].  The  load  on  the  indenter  is  then  obtained  by 
integration  as  indicated  in  Eq.  (21).  The  calculations  of  the  load  and  the  function  K{u)  for  the  given 
function  g(\)  (see  Appendix  A,  Eqs.  (A2)  and  (A  10))  are  also  carried  out  numerically  by  using  a 
Cray  supercomputer. 

For  the  conical  and  the  spherical  indenters,  Eq.  (19)  is  first  solved  by  setting  the  parameter  7 
equal  to  1  and  then  iterating  until  a  proper  value  for  7  is  obtained  such  that  //(I)  =  0.  Then  the 
correct  7  values  are  used  to  solve  the  function  H(t)  and  the  load  p . 

The  papers  of  Lebedev  and  Ufliand  [8],  Dhaliwal  [9],  and  Chen  and  Engel  [12]  contain  some 
numerical  results  and  thus  provide  an  excellent  opportunity  for  comparison  with  the  results  obtained 
here.  The  kernel  K(u),  as  defined  by  Eq.  (20),  is  computed  first  for  both  an  elastic  layer  in  smooth 
contact  with  a  rigid  substrate  and  an  elastic  layer  perfectly  bonded  to  an  elastic  substrate.  The  results 
obtained  are  the  same  as  tho.se  in  Refs.  8  and  9  (Tables  B1  and  B2  in  Appendix  B).  The  computer 
program  used  to  solve  the  Fredholm  integral  equation  (Eq.  19)  was  further  checked  by  computing 
H(t)  values  for  a  flat-ended  cylindrical  indenter  pressing  on  an  elastic  layer  in  smooth  contact  with  a 
rigid  substrate  and  comparing  the  results  with  those  given  by  Lebedev  and  Ufliand  [8].  The  excellent 
agreement  (Table  B3)  provided  confidence  in  the  correctness  of  the  present  analysis. 
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When  the  layer  and  the  substrate  are  perfectly  bonded  and  the  indenter  is  a  flat-ended  cylinder, 
the  H(j)  values  calculated  by  using  the  same  set  of  ^(u)  values  (Table  B2)  show  substantial  disagree¬ 
ment  between  our  results  and  Dhaliwal’s  published  results  (Table  B4),  as  noted,  for  example,  for  the 
case  when  hla  —  0.25  and  =  0  where  values  for  H{j)  differ  by  approximately  an  order  of 
magnitude.  However,  the  results  of  Chen  and  Engel  [12]  for  the  spherical  indenter  when  a  «  R 
were  reproduced  (Table  B5).  The  details  of  these  comparisons  are  given  in  Appendix  B. 

In  the  present  calculation,  the  number  of  terms  N  in  the  Chebyshev  series  that  approximate  the 
solution  H{j)  is  taken  to  be  5  when  hla  >  0.5.  For  values  of  h la  smaller  than  0.5,  the  number  of 
terms  should  be  increased  to  get  accurate  values  for  //(t).  For  example,  in  our  calculation,  N  was 
chosen  to  be  25  when  hla  =0.1. 


Figures  4  to  7  give  the  relationships  between  the  parameter  7  (Eq.  23)  and  the  layer  thickness  h 
for  different  indenter  shapes,  layer  and  substrate  elastic  properties,  and  bonding  conditions.  The 
parameter  7  is  the  radius  of  contact  a  normalized  with  respect  to  the  Hertz  contact  radius  an  for  a 
homogeneous  half-space  of  shear  modulus  /tj  and  Poisson  ratio  v^.  For  conical  and  spherical 
indenters  with  smooth  contact  between  layer  and  substrate,  the  results  for  t;  =  0,  0.2,  1,2,  and  10 
(where  =  (1  -  j'2)Mi/(1  ~  i'i)M2)  plotted  in  Figs.  4  and  5,  respectively.  Figures  6  and  7  show 

the  relationship  between  7  and  h  la  for  a  layer  perfectly  bonded  to  the  substrate  when 
^  =  0,  0.2,  1,  2,  and  10  (where /3  =  and  yj  =  t>2  ~  conical  and  spherical 

indenters,  respectively.  These  results  show  that  7  >  1  when  the  substrate  is  stiffer  than  the  layer 
(ri  >  1  or  0  >  1),  and  7  <  1  when  the  layer  is  stiffer  than  the  substrate  (j;  <  1  or  /8  <  1).  For  a 
given  depth  of  penetration,  the  radius  of  contact  between  the  indenter  and  the  composite  medium 
increases  with  increasing  stiffness  of  substrate  materials.  For  a  layer/substrate  stiffness  ratio  (j;  or  /3) 
smaller  than  1,  the  radius  of  contact  increases  from  the  value  a^  to  a  certain  maximum  value  and 
then  decreases  asymptotically  back  to  the  value  a^  as  the  layer  thickness  increases  from  0  to  00.  The 
layer  thickness  at  which  the  maximum  radius  of  contact  occurs  is  always  less  than  the  radius  of  con¬ 
tact  itself,  i.e.,  hla  <  \.  On  the  other  hand,  when  the  substrate  is  less  stiff  than  the  layer,  the 
radius  of  contact  decreases  from  its  homogeneous  half-space  value  a^  to  a  minimum  value  and  then 
increases  toward  a^  asymptotically  as  the  layer  thickness  h  increases  from  0  to  00.  The  minimum 
value  occurs  in  the  range  where  1  <  hla  <2.  Note  that  when  the  layer  and  the  substrate  are  iden¬ 
tical  and  there  is  no  friction  force  at  the  interface,  the  radius  of  contact  is  not  always  the  same  as  the 
homogeneous  half-space,  but  it  increases  first,  then  decreases  to  a  minimum  value  and  increases  again 
toward  the  homogeneous  value  a^- 

For  the  purpose  of  setting  up  a  guideline  for  choosing  the  appropriate  layer  thickness  and  sub¬ 
strate  properties  to  determine  the  elastic  constants  of  thin  films  by  the  indentation  test,  it  is  convenient 
to  define  a  nondimensional  quantity  p  Ip  \,  which  is  the  radio  between  p ,  the  load  needed  to  penetrate 
to  a  depth  d  into  the  composite  (layered  half-space)  and  p^  the  load  needed  to  penetrate  the  same 
depth  d  into  the  homogeneous  half-space  consisting  of  the  top  layer  material.  This  ratio  p  Ip  \  is  then 
plotted  against  h  la .  For  different  shapes  of  indenter  p  Ip ,  are  given  as  follows: 


and 


^  =  j‘//(r)rfr 
P\ 

i  H{T)dT 

p,  2  ■’o 


(cylindrical  indenter), 
(conical  indenter), 

(spherical  indenter). 


(42) 
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Figures  8  through  13  show  the  plp\  vs  hla  curves  for  the  smooth  contact  condition 
(tj  =  0,  0.2,  0.5,  1,  2,  and  10)  and  for  perfectly  bonded  condition  (/3  =  0,  0.2,  0.5,  1,  2,  and 
10,  I'l  =  >>2  =  0.3)  for  the  flat-ended  cylindrical  indenter,  the  conical  indenter,  and  the  spherical 
indenter,  respectively.  When  the  layer  is  perfectly  bonded  to,  and  softer  than  the  substrate,  the  load 
p  on  the  composite  decreases  more  or  less  exponentially  from  p^  (the  load  for  the  homogeneous  half¬ 
space  with  elastic  constants  and  to  pj  (the  load  for  the  homogeneous  half-space  with  elastic 
constants  and  i^i)  as  the  thickness  of  the  layer  increases  from  0  to  oo.  On  the  other  hand,  when 
the  layer  is  stiffer  than  the  substrate,  p  increases  asymptotically  from  P2  to  p^.  When  the  layer  is  in 
smooth  contact  with  and  stiffer  than  the  substrate,  p  decreases  from  P2  to  a  minimum  value,  then 
increases  toward  p  i  as  the  layer  thickness  varies  from  0  to  cx .  When  the  layer  is  softer  than  the  sub¬ 
strate,  p  does  not  decrease  asymptotically  from  P2  to  pi;  for  P2/M1  less  than  certain  value  (for  exam¬ 
ple  0.63  for  the  flat-ended  cylindrical  indenter  as  shown  in  Fig.  14),  the  p  values  decrease  from  p^  to 
a  minimum  value  less  than  pj,  then  increase  and  approach  pi  as  the  thickness  increases. 

When  the  indenter  is  a  flat-ended  cylinder  of  radius  a ,  the  radius  of  the  contact  area  is  always 
equal  to  a  regardless  of  the  composite  medium  material  and  layer  thickness.  This  means  the  parame¬ 
ter  7  is  always  equal  to  1.  When  the  indenter  is  conical  or  spherical,  the  radius  of  the  contact  area  is 
no  longer  a  constant,  but  is  a  function  of  the  penetration  depth  d,  the  layer  thickness  h,  the  elastic 
constants  of  the  layer  and  the  substrate,  and  the  shape  of  the  indenter.  Therefore,  the  scale  of  the 
absissa  h/a,  in  Figs.  10  to  13,  are  not  the  same  for  all  the  curves.  To  compare  the  results  for  dif¬ 
ferent  conditions,  it  is  necessary  to  express  h  in  units  of  the  constant  value  instead  of  the  variable 
a.  Figures  15  to  18  show  thep/pi  vs  h/au  curves  that  are  converted  from  Figs.  9,  10,  12,  and  13 
by  using  the  corresponding  7  values  from  Figs.  4  to  7,  respectively. 


h/a 

Fig.  4  —  Variation  of  the  normalized  radius  of  the  contact  area  between  the  conical  indenter 
and  the  smooth  contact  composite  medium  7,  with  layer  thickness  h/a,  as  a  function  of 
layer/substrate  stiffness  ratio  7; 


0.0  2.0  4.0  6.0  8.0  10.0  12.0  14.0  16.0  18.0  20.0  22.0  24.0  26.0 

h/a 

Fig.  6  —  Variation  of  the  nonnalized  radius  of  the  contact  area  between  the  conical  indenter 
and  the  perfectly  bonded  composite  medium  y,  with  layer  thickness  h/a,  as  a  function  of 
layer/substrate  stiffness  ratio  0  (v^  =  k2  =  0.3) 
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S'  I  ~  normalized  radius  of  the  contact  area  between  the  spherical  indenter 

thickness  h/a  as  a  function  of 

layer/substrate  stiffness  ratio  0  (i>i  =  02=  0.3) 


f'®'  normalized  load  p/p,  for  the  flat-ended  cylindrical  indenter  with 

ayer  thickness  h  la  for  a  smooth  contact  composite  medium  as  a  function  of  layer/substrate 
stittness  ratio,  jj 
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h/a 


Fig.  9  —  Variation  of  the  normalized  load  p  Ip  \  for  the  conical  indenter  with  layer  thickness 
h  la  for  a  smooth  contact  composite  medium  as  a  function  of  layer/substrate  stiffness  ratio  rj 


h/a 


Fig.  10  —  Variation  of  the  normalized  load  pipi  for  the  spherical  indenter  with  layer  thickness 
h  la  for  a  smooth  contact  composite  medium  as  a  function  of  layer/substrate  stiffness  ratio  i; 
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0.0  2.0  4.0  6.0  8.0  10.0  12.0  14.0  16.0  18.0  20.0  22.0  24.0  26.0 

h/a 

Fig.  U  —  Variation  of  the  normalized  load  plp\  for  the  flat-ended  cylindrical  indenter  with 
layer  thickness  h  la  for  a  perfectly  bonded  composite  medium  as  a  function  of  layer/substrate 
stiffness  ratio  (vi  =  Vi  =  0.3) 


o  0.0 
<»  0.2 
*  0.5 
c  1.0 
8  2.0 
X  10.0 


0.0  2.0  4.0  6.0  8.0  10.0  12.0  14.0  16.0  18.0  20.0  22.0  24.0  26.0 

h/a 

Fig.  12  —  Variation  of  the  normalized  load  p  /p\  for  the  conical  indenter  with  layer  thickness 
h  /a  for  a  perfectly  bonded  composite  medium  as  a  function  of  layerfsubstrate  stiffness  ratio 
(3  (w,  =  i>2  =  0.3) 
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h/a 


Fig.  13  —  Variation  of  the  normalized  load  pJp\  for  the  spherical  indenter  with  layer  thickness 
h  la  for  a  perfectly  bonded  composite  medium  as  a  function  of  layer/substrate  stiffness  ratio 
0  (ki  =  v2  =  0.3) 


Fig.  14  —  Variation  of  the  normalized  load  p/pi  for  the  flat-ended  cylindrical  indenter  with 
layer  thickness  h/a  for  a  smooth  contact  composite  medium  as  a  function  of  layer/substrate 
stiffness  ratio  r/ 


Fig,  16  —  Variation  of  the  normalized  load  p  /pi  for  the  spherical  indenter  with  layer  thickness 
h  fan  for  a  smooth  contact  composite  medium  as  a  function  of  layer/substrate  stiffness  ratio  ij 
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In  comparing  these  results,  it  is  noted  that  the  curves  are  quite  similar  for  all  three  types  of 
indenters  despite  the  large  differences  in  values  of  p  and  P|  for  the  different  indenters.  This  similar¬ 
ity  allows  us  to  establish  guidelines  for  choosing  a  priori  the  approximate  film  thickness  and  substrate 
elastic  properties  for  determining  the  elastic  constants  of  the  film  within  a  given  degree  of  accuracy. 
Figures  19  and  20  show,  for  the  perfectly  bonded  and  smooth  contact  conditions  respectively,  the  film 
thickness  needed  so  that  the  absolute  values  of  (p  -  P\)/p\  are  equal  to  2%,  5%,  and  10%.  Note 
that  the  absolute  value  of  the  ratio  (p  —  Pi)/pi  has  the  same  value  as  the  ratio  of  the  impression 
modulus  (fe  ~  where  is  the  experimentally  measured  impression  modulus  obtained 

according  to  Eqs.  (30),  (34),  or  (41). 

To  determine  the  elastic  constants  of  a  perfectly  bonded  layer  by  the  indentation  test,  as  one 
may  expect,  the  most  proper  substrate  is  the  one  that  has  the  same  elastic  constants  as  the  layer  itself, 
i.e.,  /3  =  1.  Figure  19  shows  that  when  the  /3  value  is  different  from  1,  the  normalized  layer  thick¬ 
ness,  h/aff,  needed  to  obtain  a  given  accuracy  varies  almost  linearly  with  /3.  The  slope  of  the  varia¬ 
tion  is  dependent  on  the  accuracy.  The  higher  the  accuracy  required  the  steeper  the  slope  is.  When 
the  layer  is  in  smooth  contact  with  the  substrate,  the  best  candidate  substrate  is  no  longer  the  one  that 
has  the  same  elastic  constants  as  that  of  the  layer.  In  this  case,  the  substrate  should  be  slightly  stiffer 
than  the  layer.  By  assuming  Pi  =  P2  =  0.3,  t)  (in  this  case  equal  to  j8)  is  found  to  range  from  0.87 
(for  10%  accuracy)  to  0.71  (for  2%  accuracy).  For  values  of  H/oh  greater  than  approximately  3,  the 
layer  normalized  thickness  h/o/f  needed  to  obtain  a  given  accuracy  also  varies  almost  linearly  with  r). 

In  preparing  thin  films,  the  film  is  neither  perfecdy  bonded  to  the  substrate,  nor  in  perfectly 
smooth  contact  with  the  substrate.  Figure  21  gives  the  upper  bound  h/Off  vs  p|/p2  (t>y  assuming 
i/j  =  i>2  =  0.3)  relations  to  ensure  one  can  measure  the  film’s  elastic  constants  on  the  composite  by 
the  indentation  test  to  a  given  degree  of  accuracy.  Figure  21  is  obtained  by  taking  the  upper  ranges 
of  the  curves  from  Figs.  19  and  20.  This  result  suggests  that  a  good  substrate  candidate  material  is 
about  25%  stiffer  than  the  film  material  regardless  of  bonding  conditions. 


Fig.  19  —  Relationships  between  the  layer  thickness  h  /oh  and  layer/substrate  stiffness  ratio  |3 
such  that  I  (p  -  P\)/Pi  I  is  equal  to  2%,  5%,  and  10% 
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Fig.  20  —  Relationships  between  the  layer  thickness  h  /a«  and  layer/substrate  stiffness  ratio  tj 
such  that  \(p  -  pi)/pi  I)  is  equal  to  2%,  5%,  and  10% 


0.0  0.5  1.0  1.5  2.0 

r/  or  P 

Fig.  21  —  Relationships  between  the  layer  thickness  h/an  and  layer/substrate  stiffness  ratio  ij 
or  such  that  |  (p  -  pt)/P\  I  is  equal  to  2%,  5%,  and  10% 
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SUMMARY 

The  elastic  solutions  for  the  stress  fields  in  a  composite  medium  consisting  of  an  elastic  layer 
either  perfectly  bonded  or  smoothly  overlaying  an  elastic  semi-infinite  substrate  have  been  obtained 
for  the  cases  when  the  composite  medium  is  indented  by  rigid  spherical,  conical,  and  flat-ended 
cylindrical  indenters.  Numerical  results  were  presented,  and  guidelines  were  suggested  for  the  proper 
choice  of  approximate  layer  thickness  and  substrate  elastic  properties  to  determine  the  elastic  constants 
of  the  layer  within  certain  predetermined  accuracy  by  using  the  indentation  test. 

The  following  comments  may  be  offered  in  conclusion: 

1 .  For  a  given  minimum  film  thickness  and  a  given  substrate  stiffness,  the  elastic  constants  of 
the  film  may  be  determined  within  a  calculated  accuracy  for  the  cases  when  the  Him  is  per¬ 
fectly  bonded  to  or  in  frictionless  contact  with  the  substrate.  Conversely,  if  a  given  accu¬ 
racy  is  desired,  the  required  substrate  stiffness  for  a  given  film  thickness  or  the  required 
minimum  film  thickness  for  a  given  substrate  stiffness  may  be  calculated. 

2.  In  presenting  the  effect  of  film  thickness  on  indenter  load,  the  film  thickness  h  should  be 
normalized  with  respect  to  the  Hertz  radius  an  of  the  contact  area  for  the  homogeneous 
half-space  instead  of  normalizing  it  with  respiect  to  the  actual  radius  a  of  contact  area  for 
the  composite.  The  reasons  for  this  are:  the  value  a  is  not  a  readily  measured  quantity 
because  of  its  dependence  on  the  film  thickness  and  the  elastic  constants  of  each  component 
of  the  composite,  whereas  can  be  obtained  simply  by  knowing  the  penetration  depth 
and  the  indenter  geometry. 

3.  When  h /a  >4,  the  load  vs  film  thickness  for  the  composite  depends  on  the  film/substrate 
stiffness  ratio  and  it  is  not  the  same  as  for  the  homogeneous  half-space  as  suggested  else¬ 
where  [12]. 

4.  The  p  /p\  vs  h  fan  relationships  for  three  different  shapes  of  indenters— conical,  spherical, 
and  flat-ended  cylindrical— are  quite  similar  despite  the  large  differences  in  values  of  p  and 
p  I  for  the  different  indenters. 

5.  It  should  be  pointed  out  that  this  analytical  approach  is  directly  applicable  to  problems  in 
other  engineering  disciplines  such  as,  for  example,  the  design  of  large  columns  on  elastic 
foundations. 
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Appendix  A 

DETERMINATION  OF  THE  FUNCTION  g(X) 


THE  PERFECTLY  BONDED  CASE 

When  the  elastic  layer  is  perfectly  bonded  to  the  elastic  half-space  (substrate),  the  boundary  con¬ 
ditions  at  the  interface  z  =  h  art  given  by  Eq.  (8).  By  substituting  Eqs.  (1)  and  (2)  into  Eq.  (8),  the 
following  simultaneous  linear  equations  of  the  functions  /4,  (X),  (i  =  1,  2,  . . .  ,  6)  can  be  obtained: 

w  cth  w  i4i  +  w  A2  +  +  dh  w  —  ^A^  =  0, 


(/:iCth  w  —  w)Ax  +  (ki  —  w  cth  w)A2  —  cth  w  ,43  —  /44 


-  0{k2  +  w’)/45  -  0A(,  =  0, 


[2  (1  -  v{)  ~  w  cth  w]  /1|  -f-  [2(1  —  I'l)  cth  w  -  w]A2  -  A^  -  cih  w  A^  (Al) 


+  [2(1  -  J'a)  +  w]  As  +  A^  =  0, 

[w  -  (1  -  2j'i)  cth  w]  /1|  -I-  [w  cth  w  —  (1  —  2>'i)]  /42  +  cth  w  /I3  -f  A^i 


+  [(1  —  2i'2)  +  w]  As  A(,  —  0, 

where  /3  =  Hi/nj  and  w  —  }Ji .  By  solving  Eqs.  (4.3),  (5),  and  (Al),  the  function  g(K)  for  the  per¬ 
fectly  bonded  condition  can  be  expressed  as: 

+  B2W  +  fl3sh^H')  e~'*'  (B4  -I-  B5W  -I-  B^  sh^w)  sh  w 

g(X)  =  1 - ; - ; - - - - ,  (A2) 

(C^  -f  C2W^  -I-  C3  sh^H')e  -t-  (C4  -(-  Csw^  +  Bf,  sh^vv)  sh  w 

where 

fl,  =  4a  ,02(1  +  b20)0, 

B2  =  -  1  -  3b20  +  (k2  -  2bl)0^  +  b2k20\ 

B3  =  202(1  +  2ki  +  2(>2^i  0  +  ^20^)  0, 

B4  =  ki  +  (-02  +  2^>2^1  +  40,02)  ^  +  (^1  +  40,02^:2  +  20^01)0^  +  ^2^2/3^- 
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fis  =-!-(!+  4&2)  /3  +  (1  -  +  *2/3^ 


Bg  =  *1  +  (1  +  2;fei)t2)  j8  +  (*2  +  402^2*1  +  16c, +  26, 6|)  +  jk|J8^ 


C,  =4c?  (1  +  62^), 


C2  =  1  +  3*2^  -  (*2  -  26|)  0^  -  b2k20\ 


C3  =  it,  —  (1  -  2ki)b20  +  (k2  +  16c, cl  +  26,62 )  iS^  +  62^:2^^, 


C4  =  4c^  +  4ci(2c2  +  a  lit 2)  /3  —  8c,C262i8^, 


C5  =  1  +  (*2  +  262)^  -  (1  -  262)^2/3^  -  *2^^ 


c,  =  1  —  C2  =  1  -  i'2,  6,  =  1  —  2»»,,  62  =  1  -  2i'2,  A:,  =  3  -  4>'|  and  ^2  =  3  -  4i'2. 

For  the  homogeneous  half-space  consisting  of  layer  material,  i.e.,  6  —  00,  Eq.  (A2)  gives  g(X) 
=  0.  Then  by  Eq.  (20)  we  have  k(u)  —  0  and  Eq.  (6)  becomes 


J  M(X)Jo(Xr)  dX  = /(r) 


XM(X)yo(Xr)  <fX  =  0 


(0  <  r<  c), 


(a  <  r  <  00), 


and  Eqs.  (17)  and  (19)  give  an  explicit  expression  for  the  function  H(r),  namely 

T 

ff(r)  =  Fo(t)  =  1  sin  0m  (A5) 

This  result  is  the  same  as  those  obtained  by  many  authors  [Al]  for  the  homogeneous  half-space.  On 
the  other  hand,  when  6—0,  i.e.,  a  homogeneous  half-space  with  elastic  constants  H2  and  ^2,  Eq. 
(A2)  gives 


where 


g(K)  =  I  -  7f, 


(1  -  ^2)1^1 

(1  -  »'i)M2 


Substituting  Eq.  (A6)  into  Eq.  (6)  gives 


A/*(XX/o(Xr)  i/X  =  (0  <  r<  a). 
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and 


where 


\M*(K)jQ(\r)  d\  =  0  (a  <  r  <  oo). 


M*(K)  =  Af(X)/tj. 


(A8) 


By  using  Eq.  (A7)  for  and  Eq.  (7)  for  / (r),  Eq.  (A8)  assumes  the  same  form  as  Eq.  (A4)  for  the 
homogeneous  half-space  with  elastic  constants  ^2  ^nd  t'2- 

THE  SMOOTH  CONTACT  CASE 

When  the  elastic  layer  is  in  smooth  (frictionless)  contact  with  the  surface  of  the  elastic  substrate, 
the  simultaneous  linear  equations  of  the  functions  A,(X),  (i  =  1,  2,  ....  6)  can  be  obtained  by  substi¬ 
tuting  Eqs.  (1)  and  (2)  into  Eq.  (9),  thus 

(^1  cth  w  —  w)  Aj  +  (ki  —  w  cth  w)A2  —  cth  w  ^3  —  ^4  -  j3(^2  =  0, 

[2  (1  —  j^i)  -  w  cth  w]  Aj  -h  [2(1  -  i»|)  cth  w  —  w]A2  —  A3  —  cth  w  ^4 


-I-  [2(1  -  V2)  +  +  ^6  =  0, 


(A9) 


[w  -  (1  -  2vy)  cth  w]Ai  -t-  [w  cth  w  -  (1  -  2pi)]A2  -H  cth  w  A3  -t-  A4  =  0, 


[(1  -  2^2)  +  WIA5  +  Ae  =  0. 


Equations  (4.3)  and  (A9)  give 


«(X)  =  1  - 


sh^w  +  1]  {w  -t-  sh  w  ch  w) 
w  -I-  sh  w  ch  w  +  T)  (sh^w  —  w^) 


For  the  homogeneous  half-space,  Eq.  (A  10)  gives 

g(K)  =0  A  -  00 

and 

g(\)  =  1-17  h  ^0, 

which  are  the  same  as  for  the  case  when  the  layer  is  perfectly  bonded  to  the  substrate. 
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Appendix  B 

COMPARISON  OF  NUMERICAL  RESULTS 


The  following  tables  list  some  of  the  numerical  results  obtained  by  the  authors  and  those  given 
by  other  authors  [B1-B3].  Tables  B1  and  B2  show  the  numerical  values  of  the  kernel  K{u)  and  func¬ 
tion  H(t),  respectively,  for  a  flat-ended  cylindrical  indenter  on  a  layer  that  is  in  smooth  contact  to  a 
rigid  substrate.  In  these  two  tables,  the  data  on  the  second  and  the  fourth  columns  are  obtained  by 
the  present  analysis  and  data  on  the  first  and  the  third  columns  are  those  listed  on  the  first  and  the 
fourth  columns  of  Tables  I  and  11  of  Lebedev  and  Ufliand  paper  [Bl].  The  agreement  among  these 
data  shows  that  we  can  reproduce  their  results  in  this  special  case. 

Table  B3  shows  our  values  and  those  of  Dhaliwal’s  [B2]  (the  first  and  the  last  columns  on  Table 
2  for  /8  =  0  and  0.25  in  Ref.  B2)  for  the  Kernel  K(u)  for  a  flat-ended  cylindrical  indenter  on  a  layer 
that  is  perfectly  bonded  to  an  elastic  half-space.  The  good  agreement  between  these  results  provided 
a  check  of  our  formulation  of  the  function  gfw)  for  the  perfectly  bonded  condition,  which  is  rather 
complicated  to  check  otherwise.  However,  by  using  the  same  set  of  values  of  the  kernel  K(u)  shown 
in  Table  B3,  we  cannot  reproduce  Dhaliwal’s  results  for  the  H{t)  values.  These  two  sets  of  data  are 
compared  in  Table  B4.  In  some  cases  (e.g.,  hla  =  0.25),  there  are  order  of  magnitude  differences. 
Chen  and  Engel  [B3]  have  pointed  out  that  they  could  not  reproduce  all  of  Dhaliwal’s  results. 


Table  Bl  —  Numerical  Value  of  the  Kernel  K{u)  for  a 
Flat-Ended  Cylindrical  Indenter  on  a  Film  in  Smooth 
Contact  with  a  Rigid  Half-Space. 


h  la 

=  2.0 

h  la 

=  0.5 

LU" 

YSR* 

LU 

YSR 

MM 

0.58379 

2.3349 

2.33515 

0.57985 

2.09873 

19 

0.56824 

1.5513 

1.55130 

0.6 

0.54955 

0.9875 

0.98751 

0.8 

0.5248 

0.52468 

0.5753 

0.57516 

1.0 

0.4948 

0.49479 

0.3279 

0.32805 

1.2 

0.4612 

0.46114 

0.1962 

0.19621 

1.4 

0.4251 

0.42507 

0.1291 

0.12908 

1.6 

0.38783 

0.0939 

0.09396 

1.8 

0.35053 

o.mi 

0.07367 

2.0 

0.31414 

0.0602 

0.06029 

‘’Results  from  Lebedev  and  Ufliand  [Bl]. 
*  Results  from  present  report. 
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Table  B2  —  Numerical  Value  of  the  Function  H(j) 
for  a  Flat-Ended  Cylindrical  Indenter  on  a  Film  in 
Smooth  Contact  with  a  Rigid  Half-Space 


7 

h  la 

=  2.0 

h  la 

=  0.5 

LU" 

TfMI 

LU 

YSR 

0.0 

4.321 

4.3303 

0.1 

4.303 

4.3111 

mm 

4.246 

4.2584 

1.523 

1.5236 

4.153 

4.1680 

0.4 

1.518 

1.5183 

4.022 

4.0382 

0.5 

1.511 

1.5116 

3.854 

3.8693 

0.6 

1.503 

1.5036 

3.650 

3.6642 

0.7 

1.494 

1.4943 

3.416 

3.4277 

0.8 

1.484 

1.4839 

3.157 

3.1671 

0.9 

1.472 

1.4723 

2.884 

2.8919 

1.0 

1.460 

1.4598 

2.608 

2.6138 

‘'Results  from  Lebedev  and  Ufliand  [Bl], 
*  Results  from  present  report. 


Table  B3  —  Numerical  Value  of  the  Kernel  K{u)  for  a  Flat-Ended 
Cylindrical  Indenter  on  a  Film  Perfectly  Bonded  to  an  Elastic  Half-Space 


Vl  = 

0.333,  P2  =  0-250 

/3 

u 

hla  ■■ 

=  3.0 

hla  = 

=  0.25 

Dha" 

YSF? 

Dha 

YSR 

mjgm 

-0.150197 

-0.150197 

-1.802363 

-1.80236 

Kb 

-0.149575 

-0.149575 

-1.028806 

-1.02881 

WSM 

-0.147729 

-0. 147729 

-0.278302 

-0.278302 

0.6 

-0.144714 

-0.144714 

-0.085151 

-0.085151 

0.8 

-0.140617 

-0.140617 

-0.044783 

-0.044783 

0.00 

1.0 

-0.135556 

-0.135556 

-0.028933 

-0.028933 

1.2 

-0.129669 

-0.129670 

-0.020056 

-0.020056 

1.4 

-0.123110 

-0.123110 

-0.014720 

-0.014720 

1.6 

-0.116038 

-0.116038 

-0.011301 

-0.011301 

1.8 

-0.108616 

-0.108616 

-0.008969 

-0.008969 

2.0 

-0.100997 

-0.100997 

-0.007297 

-0.007298 

■SI 

-0.106266 

-0.106559 

-1.275187 

-1.27871 

■B 

-0.105844 

-0.106137 

-0.746888 

-0.749318 

wsm 

-0.104594 

-0.104885 

-0.213979 

-0.214720 

0.6 

-0.102551 

-0.102838 

-0.062722 

-0.062538 

0.8 

-0.099773 

-0.100055 

-0.029314 

-0.028821 

0.25 

1.0 

-0.096338 

-0.096613 

-0.018148 

-0.017619 

1.2 

-0.092337 

-0.092604 

-0.012580 

-0.012105 

1.4 

-0.087871 

-0.088129 

-0.009338 

-0.008939 

1.6 

-0.083047 

-0.083296 

-0.007270 

-0.006942 

1.8 

-0.077073 

-0.078211 

-0.005849 

-0.005581 

2.0 

-0.072751 

-0.072978 

-0.004818 

-0.004597 

"Results  from  Dhaliwal  [B2]. 
*  Results  from  present  report. 
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Table  B4  —  Numerical  Value  of  the  Function  H{t)  for  a 
Flat-Ended  Cylindrical  Indenter  on  a  Film  Perfectly  Bonded 
to  a  Elastic  Half-Space. 


i>i  =  0.333,  ^2  =  0.250 

h  la  - 

=  3.0 

hla  = 

0.25 

Dha“ 

YSR'’ 

Dha 

YSR 

0.0 

1.461578 

1.40430 

119.6635 

10.7832 

0.1 

1.461143 

1.40392 

114.7609 

10.7302 

0.2 

1.459840 

1.40277 

103.0672 

10.5697 

0.3 

1.457678 

1.40087 

89.7998 

10.2972 

0.4 

1.454671 

1.39823 

77.4355 

9.90395 

0.00  0.5 

1.450839 

1.39486 

67.27787 

9.37582 

0.6 

1.446208 

1.39079 

57.93688 

8.69071 

0.7 

1.440806 

1.38604 

49.18944 

7.81816 

0.8 

1.434670 

1.38065 

40.61215 

6.73197 

0.9 

1.427838 

1.37464 

31.97879 

5.45602 

1.0 

1.420351 

1.36806 

23.50496 

4.12980 

0.0 

1.289450 

1.25805 

4.120174 

3.03548 

0.1 

1.289188 

1.25781 

A  .025161 

3.03001 

0.2 

1.288405 

1.25711 

3.806495 

3.01305 

0.3 

1.287105 

1.25595 

3.570212 

2.98279 

0.4 

1,285298 

1.25434 

3.368143 

2.93574 

0.25  0.5 

1.282993 

1.25228 

3.197656 

2.86600 

0.6 

1.280207 

1.24978 

3.036678 

2.76393 

0.7 

1.276957 

1.24688 

2.860609 

2.61524 

0.8 

1.273263 

1.24358 

2.644312 

2.40384 

0.9 

1.269148 

1.23990 

2.368507 

2.12545 

1.0 

1.264636 

1.23586 

2.039119 

1.80981 

“Results  from  Dhaliwal  [B2]. 
^  Results  from  present  report. 


Table  B5  gives  both  Chen  and  Engel’s  results  (Table  II  in  Ref.  B3)  and  our  results  for  the 
values  p*  and  d*  for  a  spherical  indenter  (a  «  R)  on  a  layer  that  is  in  smooth  contact  with  an  elas¬ 
tic  half-space  {b  =  0. 1,  =  j/2  =  0.3).  The  agreement  between  these  two  sets  of  results  in  excel¬ 

lent  only  at  large  h  la  values.  The  reason  for  the  discrepancy  at  low  h  la  values  is  believed  to  be 
due  to  the  approximation  method  used  by  Chen  and  Engel.  In  their  analysis,  they  introduce  a  “per¬ 
turbing  term”  in  the  expression  of  the  contact  pressure.  Our  comment  is  that  when  the  elastic  prop¬ 
erties  of  the  substrate  are  quite  different  from  that  of  the  layer,  or  the  layer  becomes  thin  enough, 
then  the  difference  between  the  contact  pressure  for  the  composite  and  those  for  the  corresponding 
homogeneous  half-space  is  large  enough  such  that  the  perturbation  approach  is  no  longer  appropriate. 


It  seems  to  us  that  the  nondimensional  quantities  p*  and  6*  defined  by  Chen  and  Engel  [B3]  are 
quite  ambiguous.  According  to  their  definitions,  the  p*  and  the  6*  can  be  written  as 


P 


♦  =  _L 

P\ 


(Bl) 
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Table  B-5  —  Numerical  Value  of  the  Parameter  p*  and  8* 
for  a  Spherical  Indenter  on  a  Film  in  Smooth  Contact  with 
and  Elastic  Half-Space 


Results  from  Chen  and  Engel  [B3]. 

*  Results  from  present  report. 

P*  =P/(PiyI) 
d*  =  a  /y/ 

and 

8*  =  (B2) 

ys 

In  these  expressions,  it  is  clearly  seen  that  p*  is  not  the  load  on  the  indenter  for  the  composite  nor¬ 
malized  with  respect  to  the  corresponding  load  for  the  homogeneous  half-space,  neither  it  is  equal  to 
Off /a  as  they  claim.  The  5*  value  presented  is  not  appropriate  either  because  they  set  a  equal  to 
one,  which  is  not  correct  since,  as  mentioned  previously,  the  radius  of  the  contact  area  is  not  a  con¬ 
stant  but  is  a  function  of  the  shape  of  the  indenter,  the  penetration  depth,  the  layer  thickness,  and  the 
elastic  constants  of  both  the  layer  and  the  substrate. 
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